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ABSTRACT 

Using the Very Large Telescope's Spectrograph for INtegral Field Observation in the 
Near-Infrared (VLT/SINFONI), we have obtained repeated AO-assisted, NIR spec- 
troscopy of the six central luminous, Wolf-Rayet (WR) stars in the core of the very 
young (~ 1 Myr), massive and dense cluster R136, in the Large Magellanic Cloud 
(LMC). We also de-archived available images that were obtained with the Hubble 
Space Telescope's Space Telescope Imaging Spectrograph (HST/STIS), and extracted 
high-quality, differential photometry of our target stars to check for any variability 
related to binary motion. 

Previous studies, relying on spatially unresolved, integrated, optical spectroscopy, 
had reported that one of these stars was likely to be a 4.377-day binary. Our study 
set out to identify the culprit and any other short-period system among our targets. 
However, none displays significant photometric variability, and only one star, BAT99- 
112 (R136c), located on the outer fringe of R136, displays a marginal variability in its 
radial velocities; we tentatively report an 8.2-day period. The binary status of BAT99- 
112 is supported by the fact that it is one of the brightest X-ray sources among all 
known WR stars in the LMC, consistent with it being a colliding- wind system. Follow- 
up observations have been proposed to confirm the orbital period of this potentially 
very massive system. 
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1 INTRODUCTION 

30 Doradus (NGC 2070) is an active star-forming, giant Hn 
region in the Large Magellanic Cloud (LMC). At the center 
of 30 Dor is located the very massive star cluster R136 (HD 
38268), which is regarded as the closest visible example of a 
"super star cluster", excluding its Galactic clone NGC 3603, 
which is not surrounded by a very massive cluster halo. 

The inner arcsecond of R136, denoted R136a, is unre- 
solvable by conventional ground-based telescopes, and it was 
suspected that R136a is in fact a single stellar object of mor e 
than 3000 M (|Cassinelli et~aHll98ll , lFeitzinger et al1ll980D . 

* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, 
under program ID 076.D-0563, and on observations made with the 
Hubble Space Telescope obtained from the ESO/ST-ECF Science 
Archive Facility 

f E-mail: o.schnurr@sheffield.ac.uk 



However, careful ground-based work (both direct imaging: 
Moffat & Seggewiss 1983; Moffat, Seggewiss fc Sha ra 1985; 
and speckle interferometry: IWeigelt fc Baierlll985l ) showed 
that R136a was not a single star. High-resolution imag- 
ing by the Hubble Space Telescope (HST) confirmed that 
R136a indeed consists of individual, hot and luminous stars 
l|Campbell et al.lll992j ). while the region immediately sur- 
rounding R136 contains many dozens of massive O stars, 
many of them of spectral type Q3, the hottest of all known 
O-type stars at the time (|Massev fc Hunterlll998l ). 

It has been theorized that high stellar densities in the 
cores of very massive proto-clusters are a prer equisite for 
the formation of very massive stars (Bate 2002). Moreover, 
there is empirical evidence that the mass of the most mas- 
sive clu ster member correlates w ith the total mass of the 
cluster l|Weidner fc Kroupall2006l ). Therefore, it can be ex- 
pected that the most massive stars known are found in the 
core regions of the most massive, and densest, unevolved 
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clust ers known. Given its core density of ~ 10 5 Mqpc 3 
(e.g. iMoffat et al]|l985h . R136a is hence a prime candidate 
for harboring extremely massive stars. 

It is a very remarkable fact that even the hottest and 
most massive O-type stars seem not to exce ed ~ 60 M Q 
l|Lamontagne et al ] Il996l : iMassev et ail l2002h . More mas- 
sive stars are invariably members of a very luminous and 
hydrogen-rich subtype of Wolf-Rayet (WR) stars of the 
nitrogen-rich sequence, the so-called WN5-7h stars. Studies 
using model-atmospheres indicate that these WN5-7h stars 
are not classical, core-helium burning objects, usually identi- 
fied with the WR phase, but rather core- hydrogen burning 
objects on the the main sequence; their WR-like appear- 
ance is due to their very high luminosities, log(L /L(T)) > 6.0; 
Ide Koter et aDl 19971 : ICrowther fc Dessartlll99Sh which drive 
fast stellar winds, whose high densities give rise to the 
emission-line spectrum. As has been confirmed by weighing 
WN5-7h stars with Keplerian orbits in eclipsing binaries, the 
very high luminosities of WN5-7h stars indeed correspond to 
very high masses: The prese nt record holder, NGC 3603-A1, 
tips the balance at 116 M l|Schnurr et al.ll2008ah . 

Following the? spectral classification of 

ICrowther fc Dessarj (|l998l ). R136a contains three hydrogen- 
rich WN5h stars, and one 03f/WN6 star which is a transi- 
tion type betwe en the hottest O f stars and the least extreme 
WN5-7h stars (| Walbornlll986l ) and thus weaker-lined than 
the WN5h stars. Slightly off-centered are a cooler, evolved 
WN9h star (R136b), and another WN5h star (R136c). From 
ground-based, spatiall y unresolved (i.e. integrated ), optical 
spectroscopy of R136a, IMoffat fe Seggewissl (|l983l ) reported 
a 4.377-day WNh binary in R136a, with a diluted radial- 
velocity amp l itude K ~ 38 kms -1 ; the subsequent study by 
IMoffat et al] (|l985h confirmed this finding. The situation in 
R136a seems thus very similar to that in NGC 3603, where 
from unresolved spectros copy of the central arcsecond, too, 
IMoffat fc Niemelal l| 19841 ) had identified a 3.77-day binary 
among the three WN6h stars (this finding was confirmed 
by Moffat et al. 1985 as well). Follow -up observations of 
NGC 3603 confirm ed these results l|Moffat et alj |2004 
ISchnurr et al.ll2~008al ). revealing a second close binary, C, in 
NGC 3603, with an orbital period P = 8.9 days. Hence a 
more detailed investigation of the six WNh stars in R136a 
was warranted, since it offered the potential to increase 
considerably the number of known, very massive binary 
systems. 

We have therefore obtained, for the first time, repeated, 
spatially resolved spectroscopy of the six WNh stars in R136, 
in order to single out the 4.377-day binary and to iden- 
tify any further short-period system among our targets. To 
search for eclipsing systems (photospheric or atmospheric) 
among the stars, additional optical HST photometry of our 
target stars was extracted from publicly available, archival 
imaging data. 

Ou r study also complements that of Sch nurr et al] 
(2008b), who had surveyed the 41 of the 47 known, late-type 
WN stars in the LMC according to BAT99 located outside 
R136a, and thereby concludes the efforts of the Montreal 
hot-star group to monitor every WR star in the Magellanic 
Clouds to establish their binary status. 

In the present paper, we describe the data acquisition 
and reduction (Section 2), and the data analysis (Section 
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Figure 1. Montage of the four 0.8" X 0.8" fields as seen on the 
sky, reconstructed from the SINFONI data cube. Target stars are 
identified by the BAT99 numbers. North is up and East is left. 



3). We discuss our findings in Section 4, and close with the 
summary and conclusion in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Near-infrared spectroscopy 

Targets and their relevant properties are listed in Table 
[l] giving both the nu mbers from the BAT99 catalogue 
( Brc vsacher et al] Il999h and the older Radcliffe numbers 
fe.g. lFeast et al.lll96of ^ 

Observations were carried out in service mode at the 
Very Large Telescope (VLT) with Unit Telescope 4 under 
Program-ID P076.D-0563, between November 13 and De- 
cember 5, 2005. The observations thus cover a time span of 
~22 days. We obtained repeated K-band (1.95 to 2.45 (jm) 
spectroscopy using the Spectrograph for I Ntegral Field Ob- 
servation in the Near-Infr ared (SINFONI) (|Eisenhauer et aD 
120031 : iBonnet et~al] |2004) with adaptive-optics (AO) correc- 
tion to obtain the highest possible spatial and spectral res- 
olution. BAT99-111 (R136b) served as AO reference star. 
The field of view was 0.8" x 0.8" with a "spaxel" scale of 
12.5 mas x 25 mas. Our six target WNh stars were observed 
with four pointings (with R136al, a2 and a5 all in one field), 
defined as telescope offsets from the AO guide star BAT99- 
111. Figure [1] shows a montage of the four fields as seen on 
the sky, reconstructed from the SINFONI data cube. Stars 
here are separated at least as well as on HST images. 

Total exposure times were 150s per star and per visit, 
each organized in 2 detector integration times (DITs). Given 
the brightness of our targets (K ~ ll-12mag) and the AO 
deployment, no dedicated sky frames were taken. Other cal- 
ibrations (dark and flat-field frames, and the telluric stan- 
dard star) were provided by the ESO baseline calibration. 

For most of the data red uction steps, ESO's pipeline was 
used (cf. lAbuter et al.ll2006l ). Standard reduction steps were 
taken. The two-dimensional spectra produced by each illu- 
minated slitlet were individually extracted using iraf, and 
combined into one wavelength-calibrated spectrum per star 
and visit. A main-sequence B-type star was used for telluric 
corrections. To remove the B star's Br7 absorption which co- 
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Table 1. Target list of the observations. For easier identification, both BAT99 and Ra dcliffe numbers are given to gether with the spectral 
types (based on optical spectra), (b—v) colors and E(b — v) were directly adopted from lCrowther fc Dessartl l|l998t ). while their Johnson- V 
magnitudes were used to estimate narrow-band ti-band magnitudes, applying (v — V) ~ 0.1 mag to take into account the contribution 
of emission lines. 



BAT99 


other 


spec 


V 


(b-v) 


E(b-v) 


M v 




name 


type 


mag 


mag 


mag 


mag 


106 


R136a3 


WN5h 


13.1 


+0.20 


0.28 


-6.6 


108 


R136al 


WN5h 


12.9 


+0.18 


0.28 


-6.8 


109 


R136a2 


WN5h 


13.1 


+0.21 


0.28 


-6.7 


110 


R136a5 


03If*/WN6-A 


14.0 


+0.21 


0.28 


-5.7 


111 


R136b 


WN9h 


13.4 






-6.5 


112 


R136c 


WN5h 


13.6 






-6.3 



incides with the Br7/Hel A2.166 /an emission blend of the 
target WNh stars, a Lorentzian was fitted to the absorption 
line and subtracted from the B star's spectrum. Residuals 
were very small, and proved to be harmless in the subsequent 
analysis. Finally, science spectra were rectified by fitting a 
low-order spline function to the stellar continuum. The final 
uniform stepwidth of the spectra was 2.45 A/pixel, resulting 
in a conservative three-pixel resolving power of ~3000, and 
a velocity dispersion of ~33 kms _1 /pixel. 

2.2 Optical photometry 

We have de-archived optical imaging of R136 that was ob- 
tained in Cycle 8 as part of the HST program GO-8217 (PI: 
Philip Massey). Using the Space Telescope Imaging Spec- 
trograph (STIS), 30 pairs of short (1.1-1. 3s) images were 
obtained at various roll angles through the long-pass filter 
(F28X50LP). The image size was 28" x 50" and the spatial 
scale was 0.05" pe r pixel. For more det ails on the data, we 
refer the reader to iMassev et alj {2002). Standard pipeline 
data reduction was carried out on the fly by the HST archive 
software, and no further data handling was done on our part. 



3 DATA ANALYSIS AND RESULTS 
3.1 Spectroscopy 

The mean spectra of the six WNh stars are shown in Figure 
[3 The systemic velocity RV sys of each star was measured 
by fitting Gaussians to the reasonably symmetric, top third 
of the Hell A2.188/im emission line (below that, the line 
quickly becomes too asymmetric for this). For the WN5h 
stars, which have the strongest emission lines, this was done 
for each individual spectrum of the time series. The standard 
deviation of the series was adopted as error, <thv> an d the 
error of the mean of the systemic velocity was calculated as 
eom = ctrv/vN, with N = 9 the number of spectra. 

For the two weaker-lined stars, BAT99-110 (03f/WN6) 
and 111 (WN9ha), the fit could only be carried out for the 
averaged spectrum since the S/N ratio was better. Also these 
two stars display obvious P Cygni profiles in Hell A2.188/im, 
hence two Gaussians (one absorption, one emission) plus a 
linear continuum were fitted to the entire line. The resulting 
errors on the fits are comparably large, in particular for the 
faint 03f/WN6 star, because its very weak lines remained 
buried in the noise even in the mean spectrum. Since only 



one spectrum was used, no error of the mean was calculated. 
Values are listed in Table [2] 

All four WN5h stars display systemic velocities that are 
redder (larger) than the syst emic velocity exp ected from the 
LMC, 280 ± 20 kms" 1 (e.g. iKim et"aHll998f ). This redshift 
is most likely of the same nature as that reported for the 
optical Hen A46 86 line in other WN and 03f/WN6 stars 
in the LMC (e. g. ISchnurr et al]|2008bl ). i.e. due to radiative- 
transfer effects (jHillierll 19891 ) . The mean systemic velocity of 
the four WN5h stars is (384+4) kms -1 (standard deviation), 
whose dispersion is comparable to their individual errors of 
the mean; these stars thus have systemic velocities which 
are statistically consistent with each other. 

For BAT99-110, the emission part of the Cyg profile for 
Hell A2.188/im surprisingly displays a significantly smaller 
redshift than the Hell A2.188/xm emission of WN5h stars. In 
fact, despite the large error, BAT99-110's systemic velocity 
is in excellent agreement with that of the LMC. Either near- 
infrared Hell emission lines of 03f/WN6 stars do not display 
the intrinsic redshift of their optical counterparts (possibly 
because of optical-depth effects, since their wind is thinner 
than that of, say, WN5h stars) or BAT99-110 is indeed trav- 
elling at relatively high differential velocity along the line of 
sight towards the observer. In this case, it could be that 
110 is not actually in the core, but rather ejected from it 
towards the observer, i.e. along the line of sight (also see 
below). Thus, from a point of view of cluster dynamics and 
cluster evolution, it would be very interesting to carry out a 
more detailed RV study of the massive stars i n the periph- 
ery o f R136a to check for ejected cases (see also lBrandl et all 
120071 for a more detailed discussion of this topic). 

To identify the variable stars in our sample, relative 
radial velocities (RVs) were measured by cross-correlation, 
limited to the region from 2.134 to 2.214 /im which com- 
prises the two strongest emission lines, Br7/ Hel A2.166/xm 
and He ll A2.188/im. The iterative approach of lSchnurr et al.l 
1)20080 ) was followed: A high-S/N mean spectrum of the 
time series acted as cross-correlation template for each star. 
Then, all spectra of the time series were shifted by their re- 
spective RVs, to generate a better mean template spectrum, 
and the cross-correlation was repeated with the new tem- 
plate. Figure [3] shows time-plots of these relative RVs for 
each star. 

Even without this iteration scheme, the resulting RV 
scatter of each individual star, <7Rv,xcorr is surprisingly small; 
final values ar e listed in Table [2] Following the approach of 
ISchnurr et all l|2008bl ) to derive a posteriori error bars, we 
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Figure 2. Montage of average spectra of our six target stars, arranged by spectral type and limited to the useful spectral region. The 
two strongest emission lines, Br7/Hel A2.166 fim and Hell A2.188 /an are indicated. For clarity, the spectra have been shifted by 0.25 
flux units, starting with the second from the bottom. 

Table 2. Systemic velocities, standard deviations <tkv, and errors of the mean (eom = ohv/v'N) of our program stars, obtained by 
fitting a single Gaussian to the Hen A2.188/im line. For BAT99-110 and 111, two Gaussians (emission and absorption) were fitted to the 
P Cygni profile to the mean spectrum; therefore, both RVs are given, but no error of the mean was calculated. The listed RV scatter 
<Tj}y xcorr was obtained from cross-correlation for each individual star. Those stars used as RV reference are flagged accordingly. 



BAT99 other Spectral i?V sys eom ^RV, xcorr 

name type (kms — 1 ) (kms" 1 ) (kms — 1 ) 



106 


R136a3 


WN5h 


380 ± 


11 


3.7 




14 


RV reference 


108 


R136al 


WN5h 


374 ± 


19 


6.3 




19 


RV reference 


109 


R136a2 


WN5h 


392 ± 


17 


5.7 




13 


RV reference 


110 


R136a5 


03f/WN6 


288 ± 


63 




em. 


33 


very weak P Cyg 








52 ± 


16 




abs. 






111 


R136b 


WN9ha 


348 ± 


36 




em. 


26 


weak P Cyg 








157 ± 


6 




abs. 






112 


R136c 


WN5h 


389 ± 


21 


7.0 




28 


variable? 



proceeded to construct a RV reference star by using the 
least variable stars in our sample, the three WN5h stars 
in R136a. For each of these stars, relative RVs obtained 
by cross-correlation were normalized to RV — to correct 
for different velocities of the cross-correlation template, and 
combined. For the resulting 33 data points in total, the over- 
all standard deviation ctrv, a ii = 14 kms -1 was obtained. 
This value was adopted as measurement error for each indi- 
vidual data point. 

To determine whether or not the RV variability of a 
given star is statistically significant, we adopt y/2 x o"Rv,aii = 
a- cu t as a rough estimate of the threshold above which a 
star can be considere d to be variable at the 99.9% level (see 
ISchnurr et alll2008bl for a more detailed discussion) . We ob- 



tain <T cut = 20 kms - from our three reference stars, and 
find that R136a5, R136b, and R136c exceed this threshold. 
However, the large scatter for the first two stars can be ex- 
plained by very weak lines and thus relatively lower S/N, 
which greatly affects the precision of the RV measurement. 
In contrast to this, R136c is a stronger-lined WN5h star, and 
well isolated, hence its large RV scatter is likely intrinsic. 

Including R136c in our RV reference would yield a some- 
what larger error, ORv.aii = 19 kms -1 , and accordingly 
(Tout = 27 kms -1 . Hence, R136c would have exceeded the 
cut-off value in any case (if just), thus regarded as variable, 
and been removed from the reference group. 

A period search was carried out on the RV data of 
R136c, constrained to the range from 2 to 44 days. The 
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Figure 3. Relative radial velocities as obtained from cross-correlation, normalized so that the average RV = kms -1 . Error bars (±<r) 
shown in the upper left of each panel are adopted from the reference stars, cr^y a n = 14 kms -1 . Zero velocity is shown by the dashed 
line. Left graphs show RVs plotted against time, right graphs show the data folded into the phase corresponding to the period of P = 
4.377 days of the putative binary (with arbitrary zero phase). BAT99 numbers are given in the upper right corner of each panel. No 
coherent RV variation can be seen in the data, with the possible exception or BAT99-112 (R136c; top left row). See text for more details. 



lower limit of this range is set by the mean sampling fre- 
quency, while the upper limit is twice the total time cov- 
ered by the observations, 22 days. However, due to the low 
level of variability, no significant (at or above 3ct) period 
can be found. Furthermore, the scarcity of data points leads 
to many empty phase bins if a phase-dispersion method is 
applied. Thus, more and better data are required to confirm 
the variability of R136c, and eventually its binary status. 



However, if a sine wave is fitted to the RV data of R136c 
(i.e. assuming zero orbital eccentricity), the (marginally) 
best fit is obtained for a period P = 8.2 days, with an am- 
plitude K = 42 kms -1 (see Figure For this solution, 
ff(o — c) =8 kms -1 , which is unrealistically small, and es- 
sentially produced by only one data point (around tp = 0.3), 
which does not lie on the curve. If this point is excluded from 
the fit, then <r(o — c) ~ 2 kms -1 , which is an indication for 
an ill-constrained fit (effectively six independent data points 
for four parameters) rather than real. While it is not impos- 
sible that SINFONI delivers such good data, it sheds some 
doubt onto the reality of the orbital solution. We therefore 
consider R136c as a marginal case for supporting a binary, 
subject to future confirmation. 




Figure 4. Orbital solution for BAT99-112 (R136c) if a sine wave 
(i.e. a circular orbit) is forced to the RV data obtained through 
cross-correlation. 



3.2 Photometry 

Crowded-field, quasi-white light photometry was extracted 
from each of the STIS frames using the daophot ii pack- 
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age w ith the addit i onal ALLSTAR and ALLFRAME subrou- 
tines (|Stetsonlll99l IStetsonlll994l ). The point-spread func- 
tion (PSF) was iteratively fitted using a Moffat function 
(f3 — 1.5) whose full-width at half- maximum (FWHM) was 
typically ~ 1.7 pixels. The fitting radius around stars was 
set to 4 pixels and the extraction radius of the aperture was 
1.5 FWHM ~ 2.5 pixels. Sky was determined from an an- 
nulus between 7 and 9 pixels from the center of the star. 
In order to separate all the stars in the crowded center of 
the cluster, profile-fitting was done iteratively. That way, we 
were able to obtain reliable photometry of BAT99-108, 109 
and 100, as well as 3 other fainter stars in the crowded field 
around BAT99-109. 

To carry out differential photometry, a single adjust- 
ment to the photo metric zero point was applied as in 
iMassev et all ((20021 ). For each star, the dispersion around 
the time-averaged magnitude is plotted against its instru- 
mental magnitudes (Figure [5} • Indicated are both our pro- 
gram stars and the binaries and variable stars reported by 
IMassev et all (|2002T i. Since the exposures were very short 
(~ 1.3 sec), our target stars, which are the brightest objects 
in the field, remain well below the non-linearity limit. 

These very short exposures times have repercussion on 
the photometric error statistics. If there were no exterior 
noise sources, photometric accuracy would be simply de- 
scribed by Poissonian statistics, and related to the stellar 
intensity I such that ai = VI. In reality, however, read-out 
noise becomes quickly very important for the photometric 
accuracy of fainter stars. The photometric precision of very 
bright stars, on the other hand, is limited by the quality 
of the flat-field frame. To obtain a reliable error statistics, 
we therefore have applied a more detailed noise model, tak- 
ing into account the following noise sources: i) The Poisson 
noise of the stellar intensity I, which is simply <ti. u) The 
read-out noise or, which for STIS is 4 e~ per pixel (gain = 
1). Note that the mean read-out level (R) = 0, since there is 
no flux associated with a read-out operation, only an error 
on the count rate. The extraction radius of the PSFs are 
~ 2.5 pixels (see above), hence PSFs cover a total area of 
(2.5) 2 7r ~ 20 pixels; hence the total read-out noise sums to 



Table 3. Mean instrumental magnitudes (through the long-pass 
filter) and standard deviation around the mean for our target 
stars as obtained from HST-STIS images using DAOPHOT II. 



4 x V20 = 18 e~. Hi) The noise of of the flat-field frame, 
which we assume to be normalized to unity prior to dividing 
the science exposures by it, so that the mean flat-field factor 
(F) = 1. 

Since the HST is in space and exposures are too short 
for a significant stray-light contribution, there is no sky 
background and its associated noise that need to be con- 
sidered. We further note that the flat-field noise is multi- 
plicative (it is proportional to the intensity recorded by a 
pixel) , whereas the read-out noise is purely additive. Hence, 
the photometric precision for bright stars is limited by the 
quality of the flat-field frame, whereas faint stars are more 
affected by read-out noise. 

Thus, from the usual relation 



m = m - 2.5 log(F x I + R), 

where mo is a magnitude zero point, error propagation for 
independent values of F, I, and R then yields the total error 
on the measured magnitude, 



BAT99 other STISmag a phot 
name (mag) 



106 


R136a3 


12.850 


0.013 


108 


R136al 


12.233 


0.011 


109 


R136a2 


12.863 


0.016 


110 


R136a5 


13.778 


0.017 


111 


R136b 


13.124 


0.015 


112 


R136c 


13.072 


0.013 
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(with e the Euler number), which reduces to 



Using ctp = 0.01 (which is consistent with the typical S/N 
of a flat-field frame), and mo = 24.9, this noise model fits 
the data very well. Both the la and 2a curves, based on 
this model, are shown in Figure [5] As can be seen, even the 
brightest stars in our sample do not yet reach the noise floor 
that is determined by the flat-field noise. 

As can clearly be seen, our target stars do not display 
any significant photometric variability; in fact the dispersion 
around their respective mean magnitudes is very small (see 
Table [3] for individual values). Since we were perfectly able 
to identify the eclipsing b inarie s and other variable stars 
reported by IMassev et al.l l|2002l ). we have to conclude that 
our target stars appear to be photometrically constant. This 
is especially relevant for BAT99-112 (R136c), which despite 
its short spectroscopic period does not show any indications 
of phase-dependent variability. 



4 DISCUSSION 

With our data, we cannot confirm the existence of the 4.377- 
day binary that wa s repo rted in optical em i ssion lines by 
iMoffat fc Seggewissl (|l983l ) and lMoffat et all l|l985l) . More- 
over, we cannot clearly identify any short-period binary 
among our six target stars; there only is one binary candi- 
date, BAT99-112 (R136c). This star was already known to 
be a very bright X-ray source. Several studies have reported 
a very high X- ray luminosity of Lx > 8.5 x 10 34 ergs s - 1 



for BAT99-112 rtPortegies Zwart et alj|2002l ; [Townslev et all 
120061 ; iGuerrero fc Chull2008l ). in fact the second-brightest 
object in the greater R136 area after BAT99-116 (Mk 42). It 
is very like ly that both sta rs are colliding - wind b inaries (cf. 
IUsovl ll992l ). Interestingly, ISchnurr et all l|2008bl ) reported 
significant RV variability for BAT99-116 as well, but were 
unable to establish a periodicity. While it is thus very likely 
that both are long-period systems, we note that short-period 
systems, i.e. with a small orbital separation, are not ruled 
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Figure 5. Photometric dispersion of the times-series of each source in the STIS field-of-view, plotted against its instrumental magnitude. 
The solid line is the la curve, the dashed line the 2a curve according to the model described in the text. Our target stars are identified 
by their respective BA T99 number. We have also identified eclipsing binaries (squares) and other variables (triangles) that were reported 
bv lMassev et al.l [120021) , using their nomenclature (R136-nnn). 



out by the high X-ray luminosities. NGC3603-C is an 8.9-day 
binary with o ne of the highest X-ray fluxes known a mong 
all WR stars l|MofTat et al.ll2002l ; ISchnurr et al.l [2008a! ) . so 
it seems that self-absorption of X-ray photons in the dense 
WN wind is not too much of a problem even in relatively 
close systems. 

Thus, the following question arises: If one or more of 
our target stars in R136 are indeed short-period binaries, 
could their non-detection be the result of low orbital incli- 
nation angles? After all, very low inclination angles could 
comfortably explain why neither clear RV nor photometric 
variability (i.e., eclipses or ellipsoidal variations) are found, 
despite any possible enhancement of short-period systems 
(see above). To investigate this in more detail, we have cal- 
culated the RV scatter orv, wn = K/y/2 of the primary 
(WN-type) component that would be expected from the 
continuous sampling of a circular binary with orbital pe- 
riods ranging from 1 to 100 days. In this case study, the 
masses of the two components have been fixed to 90 Mq 
for the WN5h primary and 30 Mq for the presumed O-type 
secondary, i.e. the system has a total mass of 120 Mq (the 
exact value is not critical), and a mass ratio q = M1/M2 = 3. 
The large q value is reasonably pessimistic, given that such 
large ratios are observed even in very massive systems: In 
many reported ca ses, the companion is too f aint to be (eas- 
ily) dete cted (e.g. | Schweickhardt et a.1.1 1 19991 ; Schnurr et al. 
2008a,b; ISchnurr et al.ll2009l ). 



Figure [5] Also indicated are, with two vertical, dotted lines, 
the position of the expected 4.4-day binary, and the longest 
period we can reasonably hope to find in our data set, twice 
the time coverage of our observations, i.e. ~44 days. The 
dash-dotted, horizontal line marks the observed RV scatter 
for BAT99-112, <tr V ,ii2 = 28 kms~\ 

From Figure [6] it becomes clear that if BAT99-112 were 
indeed a short-period binary, it would have to be seen under 
a very low inclination angle, i < 15°. However, in a sample 
with randomly distributed inclination angles, there would 
be 



f 90 • ■ j ■ 
J 15 o smidi 
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cos 15° 
— cos 15° 



28 



The expected ctrv,wn value for different values for the 
inclination angle, i — 15° to 90° in steps of 15°, are shown in 



times more binaries between 15° and 90° than between 0° 
and 15°; thus, it is very unlikely that there is a significant 
number of such low-inclin ation systems in our s ample . 

We not e here that IMoffat fc Seggewissl (Il983l ) and 
iMoffat et all (|l985l ) had reported a RV amplitude of K ~ 37 
kms - for R136a, which seems surprisi ngly consistent with 
Q"RV.ii 2 x t/2 = 40 kms -1 . H owever, IMoffat fc Seggewissl 
l|l983h and lMoffat et al l l|l985l ) only had integrated, i.e. spa- 
tially unresolved spectroscopy, and they reported the binary 
to be member of R136a, i.e. one of the three WN5h stars. 
(The fourth star, BAT99-110, is 03f/WN6 and hence too 
faint both in magnitude and emission-line strength to sig- 
nificantly have contributed to the integrated signal.) Our 
SINFONI observations, however, have individually resolved 
the stars; hence one would expect an RV amplitude at least 
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Figure 6. RV scatter <trv = of the WN component in 

a fictuous, continuously sampled, circular binary, with masses 
fixed at 90 M Q for the WN5h primary and 30 Mq for the O- 
type secondary, as a function of the orbital period of the sys- 
tem. Curves show different values of the orbital inclination an- 
gle, i = 90, 75, 60, 45, 30, 15°, from top. The dot-dashed, horizon- 
tal line indicates the measured scatter for R136c (BAT99-112), 
(try = 28 kms -1 . The two dotted, vertical lines indicate the 
period of the putative binary, 4.4 days, and the upper limit of 
periods we can reasonably hope to detect in our 22-day long cam- 
paign, 44 days. 

three times higher, i.e. K ~ 120 kms" 1 . Even if one allows 
for the fact that due to non-continuous sampling, the mea- 
sured <trv could be somewhat smaller, such large RV scatter 
is simply not seen among the WN5h stars in R136a. 

Hence, we have to conclude that the large scatter in 
Moffat & Seggewiss's (1983) and Moffat et al.'s (1985) data 
was not due to the binary motion of one of the target stars, 
but of different origin. A possible explanation could be that 
during the observations, the slit was not located at exactly 
the same position each time, e.g. due to difficulties intro- 
duced by variable seeing conditions, and that thus the light 
took different optical paths through the instrument. This 
could lead to an increased RV scatter, while well-isolated 
single stars are much less affected by this problem (see e.g. 
ISchnurr et afll2008bl ). It thus seems that with the exception 
of the 8.2-day (but possibly longer-period) binary candidate 
BAT99-112, there is no binary among our six program stars. 

For BAT 99-111 (R136b), this resul t does not come 
unexpectedly: ICrowther fc Dessartl (| 19981 ) re-classified this 
star as WN9ha, and there seems to be a dichotomy be- 
tween WN6,7 and WN8,9 stars when it comes to their bi- 
nary status: No WN8,9 star is known to reside in a short 
(P < 200 d ays) WR+O binary s ystem (Milky Way: iMoffatJ 
1 19891 : LMC: ISchnurr et alj|2008bh . Also, if BAT99-111 is in- 
deed a WN9ha star, then it is a descendent of an 08Iaf star 
|Crowther fe Bohannan| [l997). i.e. less massive but slightly 
more evolved than the 03f/WN6 and WN5h stars in R136. 
This could be an indication that BAT99-111 is not a member 
of R136a's relatively unevolved stellar generation, but rather 
from an older population of massive stars outside R136, and 
simply located in the line of sight. 

If, on the othe r hand, BAT99-111 is ind eed an 04Iaf+ 
star, as classified by Ma ssev fc Hunte"rl l|l998h , then both the 
spectral type and the fainter visual magnitude are consistent 



with BAT99-111 being a less extreme version of the brighter 
03f/WN6 stars in and around R136, and a member of R136a 
after all. Unfortunately, our data do not allow us to settle 
this classification issue. 

For the other five, less evolved objects in our sample, 
in particular for the three very crowded WN5h stars in the 
core of R136a, the null-result is somewhat surprising. We 
had expected to find at least one short-period binary in 
R136a, similar to the situa tion in its almost perfe ct, Galac- 
tic twin NGC3603, where ISchnurr et"all (|2008aT ) reported 
that of the three central WNh s t ars, tw o are binaries with 
P < 10 days. Also. lSchnurr et al.l l|2008bl ) report that among 
the nine 03f/WN6 and WN5-7h stars in 30 Dor, but outside 
R136, only two (~20%) are confirmed spectroscopic binaries. 
This value (at least) should have been recovered in the core 
of R136 as well, since in the present study, we are using data 
of almost identical quality. 

Considering the small number of studied stars, how- 
ever, our non-detection of systems is fully consistent with a 
frequency of ~20% of binaries w ith P < 10 days. In fact, 
combining our results with that of ISchnurr et~all ((2008bh . we 
can confirm that the short-period binary frequency among 
very massive stars around R136 within 30 Dor is not very 
high: Out of the total of fifteen 03f/WN6 and WN5-7h stars 
now studied, only two are confirmed binaries with P < 10 
days, i.e. 2/15 = 13%. If we include the candidate BAT99- 
112 (R136c) as a positive detection, we obtain 3/15 = 20% 
binary frequency. If the distribution of orbital periods were 
flat in logP (or proportional to 1/P), something which is 
usually assumed, we would find 13% (20%) binaries in each 
period bin [1;10], [10;1000], and [100;1000], i.e. the total fre- 
quency of binaries with orbital periods shorter than 1000 
days thus would be 40% (60%). 

Very recently, iBosch et alj (|2009t ) have reported on 
the results of a spectroscopic monitoring campaign of 
absorption- line OB stars in the 30 Dor region; for 52 stars, 
six to seven spectra per star were obtained over a time 
span of ~500 days. While this small number of spectra does 
not al low one to establish an orbital solution, IBosch et al.l 
(2009) were able to identify binaries from their significantly 
large RV scatter compared to single stars and binaries with 
much longer periods; among the 52 stars monitored, 25 stars 
showed RV variabi lity. If we suppose that with their method, 
IBosch et all {2009) were able to detect binaries with periods 
up to ~1000 days (i.e. twice the covered time span), the bi- 
nary frequency is thus ~ 48%. Taking into account all sta- 
tistical errors, our null-result for P < 10 days reported here 
is perfectly consisten t with this v a lue. R emarkably, Monte 
Carlo simulations by IBosch et al.l l|2009h report that their 
observational results are consistent with a binary frequency 
of 100%. Although our observations were only designed to 
identify the reported 4.4-day culprit, and we have not yet 
enough data points to identify systems much beyond P ~ 40 
days, we can rule out such a high binary frequency for our 
stars with strong winds, because only BAT99-112 shows an 
X-ray excess. 



5 SUMMARY AND CONCLUSION 

We have obtained for the first time, spatially-resolved, re- 
peated, low-spectral resolution, near-infrared spectroscopy 
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of the most luminous stars in the core of R136 in an at- 
tempt to identify th e 4.377 -day binary that w a s repo rted by 
iMoffat fe Seggewiss! (jl983r i and iMoffat et all (|l985l ). Addi- 
tional archival HST-ST1S imaging was used to extract pho- 
tometry of our target stars. 

For none of the studied stars could significant photo- 
metric variability be found. Furthermore, we cannot confirm 
the presence of a binary system with the reported 4.4-day 
period, nor do we identify any other short-period (P < 44 
days) system among the most luminous stars in R136. One 
star, however, BAT99-112 (R136c), shows small, marginally 
significant RV scatter. A forced sine-fit to the data yields a 
best period P — 8.2 days. While it is not entirely impos- 
sible, it is not very likely to see a binary under very low 
inclination angle (i < 15°), which would be required to rec- 
oncile BAT99-112's 8. 2d binary nature with the observed, 
low RV scatter. Thus, it is more probable that BAT99-112 
is a longer-period system that the limited time coverage of 
our observations could not detect, as is suggested by the fact 
that the star is one of the brightest X-ray sources among all 
known WR stars in the LMC. 

To settle this issue, long-term monitoring of BAT99- 
112 will be required, but the possibility that it is indeed bi- 
nary clearly harbors the potential to weigh one of the most 
massive stars known, and to obtain a clearer picture of the 
very upper main sequence in high-density, high-mass envi- 
ronments like R136. 

The question as to why R136 appears to have a dearth 
in its short-period binary content, compared to its Galac- 
tic counterpart NGC 3603, remains to be explained, if not 
simply due to small numbers. A key difference, though, is 
the fact that NGC 3603 is not surrounded by a massive stel- 
lar halo as is R136 (i.e. 30 Dor). Another difference is the 
factor-of-two lower metallicity in R136. How these or other 
as-yet unknown factors play out remains a mystery. In fact, 
ultimately, long-term monitoring of all the luminous stars in 
R136 is necessary to check for long-period systems. 
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